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Abstract 

We discuss and clarify a simple and accurate interpolation scheme for the spin- 
resolved electron static structure factor (and corresponding pair correlation func- 
tion) of the 3D unpolarized homogeneous electron gas which, along with some an- 
alytic properties of the spin-resolved pair-correlation functions, we have just pub- 
lished [1]. We compare our results with the very recent spin-resolved scheme by 
Schmidt et aJ. [2,3], and focus our attention on the spin-resolved correlation ener- 
gies and the high-density limit of the correlation functions. 
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1 Introduction 

The homogeneous electron gas is a model solid whose positive ionic charges 
are smeared throughout the whole crystal volume to yield a shapeless, uniform 
positive background (whence the nickname of jellium). The model, by ignoring 
the ionic lattice which makes real materials different from one another, allows 
the theorists to concentrate on key aspects of the electron-electron interac- 
tion. It thus represents an obvious limit for the inhomogeneous electron gas, 
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and, through the Density Functional Theory (DFT) [4] , its Local Density Ap- 
proximation (LDA) and other semi-local [5] and nonlocal [6] approximations, 
it links to a popular and very successful description of real materials. 

The pair-distribution functions g aiU2 (ri, r 2 ) describe the spatial correlations 
between electron pairs of prescribed spin orientations: the expected number 
of spin-cx 2 electrons in the volume dV at r 2 , when another electron of spin <Ji 
is at r 1; is equal to rfA^(r 2 cr2|r 1 cTi) = n a2 (r 2 ) g aia2 (r 1 , r 2 ) dV, where n a {v) is 
the density of spin-a electrons. In the spin-unpolarized jellium, the electronic 
spin density nj(r) = nj(r) = n/2 = (87rrf/3) _1 is uniform in space (i.e. 
independent of r), so g ai a 2 ( r i) r 2) only depends on the distance between the 
two electrons r = |ri — r 2 |. Hartree atomic units are used throughout this 
work. The static structure factor S(q) is an "experimental" quantity which 
gives a measure of the instantaneous density correlations in the system, and 
is directly related to the Fourier transform of the pair-correlation function. 
For an unpolarized homogeneous electron gas, after introducing the Fermi 
wavevector qp = (3n 2 n) 1 / 3 , the scaled variables p = qpr and k = q/qF 
are often convenient. With these variables the spin-resolved static structure 
factors are written as 

Sa x o 2 {k]r s ) = 5 aija2 + f dpp 2 [g ai(72 {p;r s ) - 1] Sm k p . (1) 

6 7T JO k p 

the total pair-distribution function is equal to g = \{g]] + gn), and the total 
static structure factor to S = + Sfj,. 

The pair-distribution functions of the uniform electron gas are a key ingre- 
dient in the construction of semi-local and nonlocal density functionals [5,6]. 
At the densities of interest for DFT calculations, the best estimate for the 
pair-correlation functions and static structure factors[f] of jellium is given by 
quantum Monte Carlo simulations (QMC) [7,8] which are available for a dis- 
crete set of interelectronic distances p (or momentum transfer k) and densities 
r s . In a recent work [1] we have presented simple functional forms for the spin- 
resolved pair-correlation functions which depend upon p and r s , are analytic 
and closed-form both in real and reciprocal space, fulfill most of the known 
properties of their exact counterparts, and contain some free paramaters that 
have been fixed by a two-dimensional (p, r s and k, r s ) fit to the new QMC 
data [8], thus yielding very accurate and reliable functions in the relevant 
density range 0.1 < r s < 10. As a byproduct, we also obtained accurate 
spin-resolved correlation energies which fulfill the exact high-density limit by 
construction [1]. In Sec. 2 we summarize the small-p (k) behavior of g uia , 2 
(S aia2 ) together with the corresponding large-/c (p) behavior of S aiU2 {g ai u 2 )i 
and we discuss and clarify some points, especially Eq. (12) of Ref. [1]. Sec. 3 is 

2 The pair-correlation functions and static structure factors are independently ex- 
tracted by quantum Monte Carlo simulations [7]. 
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devoted to a brief comparison of our results with the new spin-resolution [3] of 
the Perdew-Wang function [2], and in Sec. 4 the corresponding spin-resolved 
correlation energies are discussed. 



2 Behavior of the spin— resolved correlation functions for small and 
large arguments 



The pair-distribution function g = g cx + g c (and correspondingly the static 
structure factor) can be divided into an exchange-only contribution g e x(p) 
(given by the Hartree-Fock approximation) and a Coulomb-correlation con- 
tribution g c (p;r s ), which, in turn, can be split into its f! an d TT parts, 
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functions for small and large arguments are summarized below. 
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Here cj p = y3/r^ is the classical plasma frequency and the constants a, 6, c 
and are not known. The small-p (and corresponding large-A;) behavior of the 
correlation functions is well known [9] from the many-body Schrodinger equa- 
tion when two electrons approach each other (cusp conditions). The small-/c 
(and corresponding large-p) behavior of the |j and f"f correlation functions 
seems, instead, to be less known [1]. It can be determined by means of the 
random phase approximation (RPA) (see e.g. Ref. [10]), which is exact in the 
k — > limit, as follows. The RPA only takes into account direct processes, i.e. 
processes that occur for both parallel- and antiparallel-spin pairs. Thus, for 
an unpolarized gas, the RPA fj and f f correlation functions are equal. This 
simple consideration, together with the well-known small-A; behavior of the 
total S, qpk 2 /2u p , tells us that, as k — > 0, the linear term 3k/4 of S cx must 
be canceled 50% by and 50% by S^. The same argument can be applied 
to the — k 3 /16 term of S cx ; in fact, no term oc k 3 must appear in the total S, 
whose long-wavelength behavior is determined by the plasmon contribution 
and by the single-pair and multipair quasiparticle-quasihole excitation con- 
tributions, proportional to k 5 and k 4 respectively [10,11]. While the leading 
—3k/8 term of S% (which corresponds to the large-p term |p -4 ) and the 
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plasmon contribution oc k 2 must hold beyond RPA, the k 3 /32 term (which 
corresponds to the large-p term |p -6 ) holds for the exact S° 1(T2 in the high 
density limit, but its validity at lower densities must be verified. Notice also 
that the k 5 term in the k — > expansion of the total static structure factor 
implies that the total g behaves like p~ s when p — > oo. The functional forms 
proposed in Ref. [1] exactly fulfill all of the above analytic constraints. 



3 Spin-resolved correlation functions 



In Ref. [1] we compared our correlation functions with the widely used Perdew- 
Wang (PW) [2] model. The PW function turned out not to be accurate in its 
spin-resolved version, mainly because it does not fulfill the exact k — > limit 
of S° . Based on our work, Schmidt et al. [3] have recently proposed a new 
spin resolution of the PW model, obtained by imposing this exact limit. Such 
revised PW function represents a considerable improvement over the original 
one, and works very well for r s < 2, both in real and reciprocal space. At lower 
densities, however, our functions [1] provide a better interpolation of the QMC 
data [8], as shown, for example, by Fig. 1, where g+« and £*, for r s = 4 ob- 
tained from our scheme [1], the revised PW model [3] and QMC simulations [8] 
are reported. We have checked that, as r s increases, the discrepancy between 
revised PW and QMC data becomes more pronounced. The reason of such 
increasing discrepancy with descreasing density can be explained as follows. 
The revised ft PW function for the unpolarized gas is built up by rescaling the 
pair-correlation function of the fully polarized gas in such a way that the ex- 
act k — > limit of is fulfilled. As a result, the parallel-spin contribution to 
correlation tends to be overestimated (and consequently the ff part is underes- 
timated). Correlations are highly dominated by antiparallel-spin interactions, 




Fig. 1. Parallel-spin contribution to the pair-correlation function (left), and an- 
tiparallel-spin contribution to the static structure factor (right) for r s = 4. Solid 
line: functions from Ref. [1]; dashed line: new spin-resolved PW model [2,3]; crosses: 
QMC data [8]. 



4 



where available (see e.g. Ref. [10]), and hence simple scaling arguments which 
connect correlations in the fully-polarized gas (where fj interactions are to- 
tally absent) and ft correlations in the unpolarized gas (where | j interactions 
are present and tend to dominate the electronic correlations) will provide less 
and less reliable results at lower and lower densities, as the role of Coulomb 
correlation with respect to exchange becomes more and more important. The 
same argument applies to the correlation energy, as we shall see in the next 
Sec. 4. Notice that the overestimate of the PW ft pair-correlation function 
makes the total (exchange + correlation) parallel-spin pair-distribution func- 
tion be slightly negative near r = for densities r s > 6. Nonetheless, the new 
scaling law proposed by Schmidt et al. [3] does much better than any previous 
one. 



4 Spin— resolved correlation energies and the high density limit 

Correlation energy. The spin-resolved contributions to the correlation energy 
are defined as 

= f| £ dr' s fdppg c aia2 (p,r' s ), (2) 

and hence e c = e?i + e%*. The corresponding exact high-density limit is recov- 
ered by applying the same argument which yielded the k — > expansion of 
S° 1(T2 . In the framework of RPA (see e.g. Ref. [10]), in fact, one obtains for the 
unpolarized gas = at any r s . Since in the r s — > limit RPA is exact, 
we have 

lim eS, = lim eS T = lnr s + 0(r°). (3) 

Beyond RPA (beyond orders lnr s ) e$, and are not equal because of higher- 
order exchange terms which mainly lower the \\ correlation energy. The spin- 
resolved correlation energies presented in Ref. [1] are the best-to-date estimate 
of e£ lCT2 , since they are obtained by integrating the corresponding QMC pair- 
correlation functions [8] interpolated by our g% a and S° models which also 
incorporate the exact behavior at small and large arguments. Moreover, our 
e£ 1<72 fulfill the high-density limit of Eq. (3). Previous estimates of were 
obtained by scaling the correlation energy of the fully polarized gas. The most 
widely-used scaling laws are the Stoll et al. [12] and the Perdew-Wang [2], 
which is now available in its revised form given by Schmidt et al. [3]. In 
Fig. 2 we compare our with these three schemes. It is apparent that the 
revised PW gives the best result, even if it does not fulfill the high-density 
limit of Eq. (3). As a consequence, the revised PW underestimates by ~ 22% 
the t| correlation energy at r s = 0. As expected from the corresponding ff 
pair-correlation function (see Fig. 1), at densities r s > 2 the revised PW 
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Fig. 2. Parallel-spin contribution to the correlation energy obtained from different 
scaling guesses compared to our interpolation scheme (solid line, Ref. [1]). Units are 
Hartrees per electron. 

overestimates by an amount which increases with r s (e.g. 3% at r s = 3, 
21% at r s = 10). Again, we see the failure of simple scaling laws which try 
to connect correlations in the fully polarized gas and f"f correlations in the 
unpolarized gas. 

High-density limit of the pair-correlation functions. It is also wortwhile to 
discuss some points about the link between the high-density limit of the pair- 
correlation function and the r s — > limit of Eq. (3). Rassolov et al. [13] have 
recently computed the r s — > limit of g c /r s , which turned out to be a well 
defined, r s -independent, function. For different reasons, neither PW [2,3], nor 
our pair-correlation function [1] fulfill this limit. As pointed out in Ref. [1], our 
simple functional forms do not reconcile the known high-density limit of the 
pair-correlation functions at zero interelectronic distance [14] with Eq. (3). In 
this respect, the functional form used by PW [2] is such that when the exact 
k — > limit of S^ ia2 is imposed to it, Eq. (3) is automatically violated. The 
PW spin-resolved version, thus, suffers from a different problem than ours. 
In summary, although major steps forward have been achieved in this area 
by Refs. [1,3], none of the existing models fulfills all the known properties at 
r s — > 0, and further improvements are needed in this limit. 



5 Conclusions 

In this work we have clarified the behavior of the spin-resolved pair-correlation 
function and static structure factor of the unpolarized uniform electron gas 
for small and large arguments. We have then compared the two best-to-date 
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models for these functions [1,3], pointing out their advantages and drawbacks, 
and we have discussed the corresponding spin-resolved correlation energies. 
We found that the functions of Ref. [1] provide a better fit to the QMC data [8], 
both in real and reciprocal space, and provide the best-to-date estimate of 
the spin-resolved correlation energies. At very high density, however, even the 
forms proposed by Refs. [1,3] are inadequate, and further work is needed. It 
should be also kept in mind that in this high-density limit the non-relativistic 
hamiltonian, on which the whole theory is based, is no longer valid, and rela- 
tivistic corrections or a fully relativistic treatment wolud be required [15]. 
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